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Using a data sample of integrated luminosity of about 33 pb~^ collected around 3.773 GeV with 
the BESII detector at the BEPC collider, the semileptonic decays D^ -^ (jye'^u,., D^ -^ (jy^.'^u^ and 
the hadronic decay D^ -^ (j>Tr'^ are studied. The upper limits of the branching fractions are set 
to be BF{D+ -^ (jte+Ve) < 2.01% and BF{D+ -^ -^m'^^'m) < 2.04% at the 90% confidence level. 
The ratio of the branching fractions for D'^ -^ (piv'^ relative to -D"*" -^ K~-k^-k^ is measured to 
be 0.057 ± 0.011 ± 0.003. In addition, the branching fraction for D"*" — > (/iTr"*" is obtained to be 
(5.2 ±1.0 ±0.4) X 10"^ 



I. INTRODUCTION 

Searching for new modes in charm decays is of great 
interest. It not only investigates possible decay mecha- 
nism and finds its contribution to the total decay width, 
but is also useful to simulate accurately cascade decays of 
bottom mesons and to eliminate backgrounds of charm 
modes in studying bottom decays. 

In this paper, the semileptonic and hadronic D decays 
in which the final state particles contain a </> meson are 
studied. Whenever a specific state or decay mode is men- 
tioned in this work, the charge-conjugate state or decay 
mode is always implied. 



II. BESII DETECTOR 

The BESII detector upgraded from the BES [H is a 
large solid-angle magnetic spectrometer described in de- 
tail elsewhere |3l. A 12-layer Vertex Chamber (VC) sur- 
rounding the beryllium beam pipe provides the trigger 
and coordinate informations. A forty-layer main drift 
chamber (MDC), located outside the VC, yields precise 
measurements of charged particle trajectories with a solid 
angle coverage of 85% of 47r; it also provides ioniza- 
tion energy loss (dE/dx) measurements used for particle 
identification. Momentum resolution of 1.78%-\/l + p^ 
{p in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering are obtained for data taken at ^/s — 3.773 
GeV. An array of 48 scintillation counters surrounding 
the MDC measures the time of flight (TOF) of charged 
particles with a resolution of about 180 ps for electrons. 
Outside the TOF is a 12 radiation length barrel shower 
counter (BSC) comprised of gas tubes interleaved with 
lead sheets. The BSC measures the energies of electrons 
and photons over 80% of the total solid angle with an 
energy resolution of cje/E ~ Q.22/\/E {E in GeV) and 
spatial resolution of a^ = 7.9 mrad and <Jz — 2.3 cm for 
electrons. A solenoidal magnet outside the BSC provides 
a 0.4 T magnetic field in the central tracking region of 
the detector. The magnet flux return is instrumented 
with three double layers of counters, that are used to 
identify muons with momentum greater than 500 MeV/c 
and cover 68% of the total solid angle. 



III. DATA ANALYSIS 

The data used for this analysis were collected around 
the center-of-mass energy of 3.773 GeV with the BE- 
SII detector operated at the Beijing Electron Positron 
Collider (BEPC). The total integrated luminosity of the 
data set is about 33 pb~^. At the center-of-mass en- 
ergy 3.773 GeV, the ^'(3770) resonance is produced in 
electron-positron (e+e~) annihilation. The ■!/'(3770) de- 
cays predominately into DD pairs. If one D meson is 
fully reconstructed, the D meson must exist in the system 



singly tagged D). Using the singly tagged D~ sample, 
the semileptonic decays D"*" -^ ^e+j/g and D^ -^ (jj^'^v^ 
are searched in the recoiling system. The hadronic candi- 



dates £)+ 



67r+ and D^ -^ K 7r+7r+ are reconstructed 



directly from the data sample of 33 pb 



A. Event selection 

Events which contain at least three charged tracks with 
good helix fits are selected. To ensure good momentum 
resolution and reliable charged particle identification, ev- 
ery charged track is required to satisfy |cos6'| < 0.85, 
where 9 is the polar angle. All tracks, save those from Kg 
decays, must originate from the interaction region by re- 
quiring that the closest approach of a charged track is less 
than 2.0 cm in the xy plane and 20 cm in the z direction. 
Pions and kaons are identified by means of the combined 
particle confidence level which is calculated with infor- 
mation from the dE/dx and TOF measurements |^|j|. 
Pion identification requires a consistency with the pion 
hypothesis at a confidence level {CL-^) greater than 0.1%. 
In order to reduce misidentification, a kaon candidate is 
required to have a larger confidence level [CLk) for a 
kaon hypothesis than that for a pion hypothesis. For 
electron or muon identification, the combined particle 
confidence level (CLg or CL^), calculated for the e or 
^ hypothesis using the dE/dx, TOF and BSC measure- 
ments, is required to be greater than 0.1%. 

The 7r° is reconstructed in the decay of 7r° -^ 77. To 
select good photons from the tt*^ decay, the energy de- 
posited in the BSC is required to be greater than 0.07 
GeV, and the electromagnetic shower is required to start 
in the first 5 readout layers. In order to reduce back- 
grounds, the angle between the photon and the nearest 
charged track is required to be greater than 22°, and the 
angle between the cluster development direction and the 
photon emission direction to be less than 37° [3, 3 • 



B. Singly tagged D sample 



The singly tagged D~ saniple used in the analysis 
was selected previously H, H, IE 111 - The singly tagged 
D~ mesons were reconstructed in the nine hadronic 
modes of D" -^ K+tt-tt-, K%it- , K^R- , K-K+n', 

°, iir+7r^7r^7r^7r+ and 
TT TT TT ' . The distributions of the fitted invariant masses 
of the mKm:(m = 0, 1, 2; n = 1, 2, 3, 4) combinations are 
shown in Fig. ^ The number of the singly tagged D~ 
mesons is 5321 ± 149 ± 160 H,|y]i where the first error is 
statistical and the second systematic. 



K°Tr-TT-n+, K°n-n°, K+n n ir 



C. Candidates for D'^ -^ 4>^'^^e and D^ -^ (j)jj.'^i'^ 
Candidates for D^ -^ 4>e'^Ve and D^ -^ (pji^v^ are 

fifAantf^A frnm fliA .Qnr-vriiT-incr -t-.rarV.Q in f.Vip eirQf.pm rpmilincr 
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Invariant Mass (GeV/c^) 

FIG. 1: The distributions of the fitted invariant masses 
of (a) K+-K--K-, (b) K%-K-, (c) K%K-, (d) K+K'tv-, 
(e) KI-K-TX--K+, (f) K°s-K--K°, (g) K+-K--K-Tx°, (h) 
A'^TT^TT^Tr^TT^, (i) TT^TT^TT^ Combinations; (j) is the fitted 
masses of the mKnn combinations for the nine modes com- 
bined together. 



cays with neutral pion, for example D~^ —> ^tt+tt'', the 
number of the isolated photons is required to equal two if 
the singly tagged D~ mode contains a 7r° meson, other- 
wise zero. The events with isolated photons those ener- 
gies are larger than 0.1 GeV excluded in the tagged D~ 
are not kept. 
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FIG. 2: The distributions of the invariant masses of (j)/! 
(solid) and ^vr"*" (dashed) for D^ — > tfju'^i/fj, and D'^ — > (/)7r^ 
from Monte Carlo events. 



against the tagged D^ . To select these candidates, it 
is required that there are three charged tracks, one of 
which is identified as an electron or a muon with the 
charge opposite to the charge of the tagged D^ meson, 
the other two are identified as a if + and a K^ mesons. 
The (j) nieson is selected by requiring that the difference 
between the invariant mass of K^K^ and the nominal 
mass should be less than 20 MeV/c^. 

In the semileptonic decays, one neutrino is unde- 
tected. A kinematic quantity Umiss = E^iss - Pmiss 
is used to obtain the information of the missing neu- 
trino, where Emiss and Pmiss are the total energy and 
the total momentum of all missing particles respectively. 
Monte Carlo study shows that the background modes 
for _D+ -^ 0e+(/z"'")r^e(u) come primarily from Z?+ -^ 



(j)li+{e+)v^(^), D+ -> (l)n+Tr° and D+ -^ <J)tt+ . The re- 
quirement |[/miss| < ScTjy^,.^^ is imposed to reduce these 
backgrounds, where crij,„i^^ is the standard deviation of 
the Umiss distribution obtained from Monte Carlo simu- 
lation. In Fig. 13 the solid and dashed histograms show 
respectively the distributions of the invariant masses of 
(j)li^ and (/)7r"'" for £)"*" — > 0//+i'^ and D~^ — > ^tt"^ from 
Monte Carlo events. The candidate mode D^ -^ 0/i+i/^ 
has a potential hadronic background of 13+ -^ ^tt"*" due 
to the misidentification of a charged pion as a muon. 
These backgrounds can be suppressed by requiring the 
invariant masses of the 0/1+ combinations to be less than 
1.8 GeV/c^. Similarly the invariant masses of the (pe'^ 
combinations are also required to be less than 1.8 GeV/c^ 
in order to remove the £)+ -^ (/jTr"*" background for the 
£)+ — > (Ae+i/« decav. To suDoress backgrounds from de- 



Fig. Eland Fig. 01 show respectively the scatter plots of 
the invariant masses of the K~^K~ combination versus 
the mKm: combination for D^ — > (pe^Vf, and D+ -^ 
(pfX^v^ before and after applying both Umiss and M^i+ 
cuts. No event for the D'^ -^ 4>e'^Ve and D^ -^ (pfX^v^^ 
decays satisfies the selection criteria. 
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FIG. 3: Scatter plot of the K^ K~ invariant masses versus 
the mKmr invariant masses for the D+ — ^ (pe'^u^ candidates: 
(a) before and (b) after applying both Umiss and M^g+ cuts, 
where the rectangle represents the combined signal region of 
(h and D+ . 
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FIG. 4: Scatter plot of the K'^K invariant masses versus 



4>fJ.'^ I'fj. candidates: 
and M^^+ cuts, 
where the rectangle represents the combined signal region of 
and 73+. 



the mKnn invariant masses for the D — 
(a) before and (b) after applying both Ur. 
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FIG. 5: The distributions of the fitted invariant masses of (a) 
K^K~ combination in the D'^ signal region, (b) K'^K~tv^ 
combination in the tj) signal region. 



D. Candidates for i3^ 



tt"*" and D 



+ 



K- 



m ' 



we 



combination 
4) mass, and 



The candidate for D^ — > 07r+ is reconstructed from 
£)+ -^ K^K^TT^ . The distribution of the fitted invariant 
mass of the K^K^ii'^ combination is shown in Fig.QJd). 
As mentioned in previous subsection, the cj) meson is se- 
lected though its decay to K^K^ . Fig. |3Ja) shows the 
fitted invariant mass spectrum of the K^ K~ pairs from 
the Z?"*" signal region for the D+ -^ K~^K~7t~^ candi- 
date events. Fitting the mass spectrum with a Gaus- 
sian function convoluted Breit-Wigner gives 64.3 ± 10.3 
(j) signal events. To select the decay £)+ 
require the invariant mass of the K~^K^ 
to be within 20 MeV/c^ of the nominal i 
the invariant mass distribution of K~^K~Tr^ combina- 
tion is shown in Fig. [S{b), a clear signal of the decay 
D+ -^ (f)TT^ is observed. Fitting the mass spectrum of 
the K^K^TT^ combination with a Gaussian function as 
the signal and a special function 0, |j| for the back- 
ground yields 55.6 ± 8.6 signal events, where the mass 
resolution is fixed at 0.0022 GeV/c^ determined from the 
Monte Carlo simulation. There may be the K'^K~ com- 
binatorial background in the observed ^tt"*" events. The 
background events are estimated to be 6.1 ± 2.5 by the 
(j) sideband \Mk+k- - 1-10| < 0.05 GeV/c^. For the 
Z3+ -^ K~7T^Tr'^ candidates, the fitted invariant mass 
spectrum is shown in Fig. QJa), and 3153.8 ± 69.2 events 
are yielded from the fit. 

In the fitted yields, there are still some background 
contaminations which shape a peak under the D~^ peak. 
These backgrounds are estimated by analyzing the Monte 
Carlo sample which is about 14 times larger than the 
data. The Monte Carlo events are generated as e+e" -^ 
DD. and both D and D mesons decav to all nossible 



modes according to the decay modes and the branching 
fractions quoted from PDG Q excluding decay modes 
under study. The number of events satisfying the se- 
lection criteria is then normalized to the data. For the 



D+ 
D+ 



decay, the dominant background modes are 



<i,^i^ 



D+ 



and D+ — > 07r+7r°, and the 



number of the background events is 3.6 ±1.1. For the 
£)+ -^ K~TT~^TT^ decay, the number of the background 
events is 35.2 ±4.2, After subtracting the numbers of the 
background events, 45.9 ± 9.0 and 3118.6 ± 69.3 signal 
events for the D~^ -^ (j>Tr~^ and D~^ — > K^tt^tt^ decays 
are remained. 



IV. RESULTS 



A. Monte Carlo efficiency 



The reconstruction efficiencies of the semileptonic de- 
cays _D+ ^ ipe^Ve, D^ — > (f)^'^v^ and the hadronic 
decays D'^ -^ 07r+, D^ -^ K^tt^tt^ are estimated by 
Monte Carlo simulation with the GEANT3 based Monte 
Carlo simulation package fl. A detailed Monte Carlo 
study shows that the efficiencies are 6^^+^^ = (2.46 ± 
0.03)%, e^^+^^ = (2.43 ±0.03)%, e^^+ = (6°47±0.05)% 
and eK-^+^+ = (24.99 ± 0.13)%. 



B. Branching fraction 

A Z)+ signal is observed neither in the D~^ — > (pe^Ve 
decav nor in the D~^ —^ 6u&v,, decav. The uDoer limit of 



the branching fraction can be calculated using 



BF < 



N. 



sg 



exN,% x{l-S)' 



TABLE II: Ratio of partial widths i'Id+->k-'w+1+) ^'^'^ com- 
Cl) parison with other experiments (Evts stands for D'^ -^ 01"^). 



where N^g is the upper limit of the signal yield given 
with the Feldman-Cousins prescription ^], which is 2.44 
for zero observed event in the absence of background for 
the confidence level of 90%, e is the detection efficiency, 
N^iag is the number of the singly tagged D~ mesons and 5 
is the systematic error. The upper limits of the branching 
fractions are summarized in the second column of Tabled 
Comparisons with those obtained by MARKIII 10] and 
BESI [ij experiments are also listed in the third and 
fourth columns of Tabic U 



Reference 






Evts 



BESII(this work) 0.057 ± 0.011 ± 0.003 46 



E687 fl2] 
WA82 [13] 
CLEO |14] 
NA14 |15] 
E691 |16] 
MARKIII fl7] 



0.058 ± 0.006 ± 0.006 - 

0.062 ± 0.017 ± 0.006 19 

0.077 ± 0.011 ± 0.005 128 

0.098 ± 0.032 ± 0.014 12 

0.071 ± 0.008 ± 0.007 84 

0.084 ±0.021 ±0.011 21 



TABLE I: Upper limits for the branching fractions (%) are 
given at the 90% confidence level. 



In addition, the branching fraction for D^ 
obtained to be 



OTT ' IS 



Mode 



BESII MARKIII BESI 



D+ 
D+ 






< 2.01 

< 2.04 



< 2.09 

< 3.72 



1.38 



The systematic error, S, includes the uncertainties from 
particle identification (1.4% for _D+ -^ (pe'^iye, 1-6% for 
£)+ — > (/)/!+ j/p), tracking efficiency (2.0% per track), pho- 
ton reconstruction (2.0%), Umiss selection (0.6%), the 
number of the singly tagged D~ mesons (3.0%) and 
Monte Carlo statistics (1.2%). These uncertainties are 
added in quadrature to obtain the total systematic errors 
to be 7.3% for D+ ^ (j)e+i^e and 7.4% for D+ -^ (j)^i+iyf,. 
The ratio of the branching fractions for £)+ -^ (pn^ rel- 



ative to D+ -^ K- 



can be obtained by 



BF{D^ 



(j)TT^ 



N, 



4>-K^ 



BF{D+ -^ K-TT+Ti^ 



^^ K 7r+7r+ 



^K 7r+7r+ 



t^TT^ 



(2) 



where N^t^+, N^- 71+11+ are the numbers of the signal 
events, and e^jr+i ^K-Tr+T!+ are the detection efficiencies 
for the tagging and normalizing modes. Inserting these 
numbers into Eq. 121), the relative branching fraction is 
obtained to be 0.057 ± 0.011 ± 0.003, where the first er- 
ror is statistical, and the second systematic. By nor- 
malizing Z?+ — > 07r+ to D+ -^ K^TT^TT^, the system- 
atic uncertainty from tracking efficiency can be canceled. 
The residual error, including both contributions of the 
tagging and normalizing modes, arises from uncertain- 
ties in particle identification (0.5% for uncanceled kaon 
track in _D+ — > (/'7r+, 0.5% for uncanceled pion track in 
D"*" -^ ii'~7r+7r+),the number of the singly tagged D^ 
meson (3.7% for D+ -^ (/)7r+, 3.3% for D+ -^ K'TT+n+), 



background uncertainty (1.4% for _D+ 
D+ ^ K 

D+ -^ (hT, 



07r+, 0.2% for 



7r+7r+) and Monte Carlo statistics (0.8% for 
-, 0.5% for D+ -^ K-TT+-K+). These uncer- 



tainties are added in quadrature to obtain the total sys- 
tematic error to be 5.3%. The relative branching fraction 
is comnared with other measurements in Table ITTl 



BF{D^ 



(5.2 ±1.0 ±0.4) X IO- 



CS) 



by using the branching fraction of the D+ — > X^tt+tt^ 
decay quoted from PDG [2|, where the first error is sta- 
tistical and the second systematic. The systematic error 
of the branching fraction for £)+ — > ^tt"*- includes the sys- 
tematic error of the ratio BF{D+ -^ (J)t:+)/BF{D+ -^ 
K~TT^TT^) ( 5.3%) and the uncertainty of the branching 
fraction for D+ -^ K-'it+tt+ (6.5%). The total system- 
atic error is obtained to be 8.4% by adding these uncor- 
related errors in quadrature. 



V. SUMMARY 

Using a data sample of integrated luminosity of 33 
P&-1 collected around 3.773 GeV with the BESII detec- 
tor at the BEPC, the semileptonic decays D'^ -^ (pe'^Ve, 
Z3+ -^ (fifi^v^ and the hadronic decay D^ -^ (/)7r+ are 
studied. The upper limits of the branching fractions are 
set to be BF{D+ ^ cj)e+Ve) < 2.01% and BF{D+ -^ 
(^/x+j/^) < 2.04% at the 90% confidence level. The ra- 
tio of the branching fractions for D"*" -^ (pTr~^ relative to 



D^ 



K- 



is measured to be 0.057± 0.011 ± 0.003. 



In addition, the branching fraction for D^ —>■ cjjtt^ is 
obtained to be (5.2 ± 1.0 ± 0.4) x lO'^. 



Acknowledgments 

The BES collaboration thanks the staff of BEPC for 
their hard efforts. This work is supported in part by 
the National Natural Science Foundation of China under 
contracts Nos. 10491300, 10225524, 10225525, 10425523, 
the Chinese Academy of Sciences under contract No. KJ 
95T-03, the 100 Talents Program of CAS under Con- 
tract Nos. U-11, U-24, U-25, the Knowledge Innova- 
tion Project of CAS under Contract Nos. U-602, U- 
34 ('IHEP'l. the National Natural Science Foundation of 



China under Contract No. 10225522 (Tsinghua Univer- 
sity). 



[1] BES CoUab, J.Z. Bai et al, Nucl. Instr. Meth. A344 [10 

(1994) 319. 
[2] BES CoUab, J.Z. Bai et al, Nucl. Instr. Meth. A458 [11 

(2001) 627. 
[3] BES CoUab, M. Ablikim et al, Phys. Lett. B597 (2004) [12 

39. 
[4] BES CoUab, M. Ablikim et al, Phys. Lett. B603, (2004) [13 

130. 
[5] BES CoUab, M. Ablikim et al, Phys. Lett. B608 (2005) [14 

24. 
[6] BES CoUab, M. Ablikim et al, Phys. Lett. B610 (2005) [15 

183. 

[7] Particle Data Group, Phys. Rev.D66 (2002) 010001. [16 

[8] BES CoUab, M. Ablikim et al, Nucl. Instr. Meth. A552 

(2005) 344. [17 

[9] G.J. Feldman and R.D. Cousins, Phys. Rev. D57 (1998) 

3873. 



Marklll CoUab, Z. Bai et al, Phys. Rev. Lett. 66 (1991) 

1011. 

BES CoUab, J.Z. Bai et al. High Energy Phys. and Nucl. 

Phys. 26 (2002) 1093. 

E687 CoUab, P.L. Frabetti et al, Phys. Lett. B351 (1995) 

591. 

WA82 CoUab, M. Adamovich et al, Phys. Lett. B305 

(1993) 177. 

CLEO CoUab, M. Daoudi et al, Phys. Rev. D45 (1992) 

3965. 

NA14 CoUab, M.P. Alvarez et al, Phys. Lett. B246 

(1990) 261. 

E691 CoUab, J.C. Anjos et al, Phys. Rev. Lett. 60 897 

(1988). 

MARKIII CoUab, R.M. BaUrusaitis et al, Phys. Rev. 

Lett. 55 150 (1985). 



